Dynamic environment coupling induce synchronized states in coupled time-delayed 

electronic circuits 
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We experimentally demonstrate the occurrence of various synchronized states in coupled piece- 
wise linear time-delayed electronic circuits using dynamic environment coupling where the environ- 
ment has its own intrinsic dynamics via feedback from the circuits. We carry out these experiments 
in two different coupling configurations, namely mutual and subsystem coupling configurations. De- 
pending upon the coupling strength and the nature of feedback, we observe a transition from nonsyn- 
chronization to complete synchronization via phase synchronization and from nonsynchronization 
to inverse synchronization via inverse-phase synchronization between the circuits in hyperchaotic 
regime. Snapshots of the time evolution, phase projection plots and localized sets of the circuits as 
observed experimentally from the oscilloscope, along with supporting numerical simulations confirm 
the existence of different synchronized states. Further, the transition to different synchronized states 
can be verified from the changes in the largest Lyapunov exponents, Correlation of Probability of 
Recurrence and Correlation Coefficient as a function of the coupling strength. We present a detailed 
linear stability analysis and obtain conditions for different synchronized states. 

PACS numbers: 05.45.Xt,05.45.Pq 



Synchronization is one of the collective phe- 
nomena most intensely studied in the chaotic dy- 
namical systems. Considerable amount of re- 
search has been focused on small assemblies of 
low dimensional coupled chaotic systems and var- 
ious types of synchronization have been demon- 
strated both numerically and experimentally. Re- 
cent studies on synchronization have also focused 
on higher-dimensional systems such as time-delay 
systems, which exhibit highly complicated hy- 
perchaotic attractors with multiple positive Lya- 
punov exponents. Different kinds of synchroniza- 
tions and their transitions have also been identi- 
fied and reported in these systems with different 
coupling configurations using numerical simula- 
tions but only limited studies are available from 
experimental point of view. In this paper, we 
experimentally investigate and report the exis- 
tence of various types of synchronizations and 
their transitions in time-delayed electronic cir- 
cuits which are coupled through a common dy- 
namic environment along with supporting nu- 
merical results. We demonstrate different syn- 
chronization transitions in two coupling configu- 
rations, namely mutual and subsystem coupling 
configurations. We have also presented a linear 
stability analysis to find the conditions for differ- 
ent synchronized states. 



INTRODUCTION 



Synchronization is a ubiquitous phenomenon often ob- 
served in coupled chaotic and hyperchaotic systems [l|, Q . 
Depending upon the strength and the nature of coupling, 
various types of synchronizations such as phase synchro- 
nization (PS), complete synchronization (CS), general- 
ized synchronization, lag and anticipatory synchroniza- 
tions, inverse phase synchronization (IPS) and inverse 
synchronization (IS) have been observed in coupled dy- 
namical systems. All these types of synchronizations 
are achieved in ensembles of dynamical systems through 
some common coupling schemes, namely linear and non- 
linear error feedback couplings [3|, mil coupling via dissim- 
ilar and/or time-delayed variables [5|-[3, inhibitory cou- 
pling [8[, coupling via dynamical relaying [9], adaptive 
coupling [lOj and also systems driven by a common noise 
[ll| etc. 

In many real world systems, synchronization can occur 
due to interaction through a common dynamic medium 
and hence the systems also evolve similar to their envi- 
ronment under the influence of the latter. Examples in- 
clude synchronization of chemical and genetic oscillators 
J12|-[l4{, synchronized behavior with self pulsating peri- 
odic and chaotic oscillations produced by an ensemble of 
cold atoms interacting with a coherent electromagnetic 
field |15l |. and in synchronizations of cells, in coupled cir- 
cadian oscillators due to common global neurotransmit- 
ter oscillation [16|. In all the cases, the coupling function 
has a dynamics modulated by the system dynamics. In 
this connection, Resmi et al [17] have recently shown and 




FIG. 1: Circuit diagram of a single time delayed feedback 
oscillator with a nonlinear device (ND) unit, a time-delay 
unit (DELAY) and a low pass first-order RqCq filter. 



verified numerically the existence of various types of syn- 
chronizations in low-dimensional chaotic systems which 
are coupled through a dynamic environment without in- 
trinsic or coupling time-delay. 

Time-delay systems are an important class of dynami- 
cal systems which occur in many real world systems and 
synchronizing such systems is having potential applica- 
tions in diverse areas of science, engineering and technol- 
ogy [18M22J due to their hyper chaotic nature. During the 
past few years, researchers are interested in synchronizing 
such time-delay systems and mostly all known synchro- 
nizations and their transitions are identified and reported 
in two coupled time-delay systems using numerical simu- 
lations [2|. In contrast, only very few studies are available 
from an experimental point of view. Experimentally gen- 
erating and synchronizing such hyperchaotic electronic 
signals (with multiple positive Lyapunov exponents) is 
very important where these signals may be used to hide 
secret messages in the area of secure communication and 
cryptography [23j. 

In this paper, we demonstrate experimentally the oc- 
currence of various types of synchronizations and their 
transitions in indirectly coupled time-delayed electronic 
circuits with intrinsic time-delay using dynamic environ- 
ment coupling. We carry out these studies in two dif- 
ferent coupling configurations, namely mutual coupling 
configuration where both the circuits and environment 
are mutually sharing their feedback, and subsystem cou- 
pling configuration where both the circuits are sharing 
their feedback with the environment, while only one of 
the circuits is receiving feedback from the environment. 
Depending upon the coupling strength and the nature of 
feedback, we observe different types of synchronization 
transitions in the coupled circuits which include tran- 
sition from nonsynchronization to CS via PS and from 
nonsynchronization to IS via IPS in both coupling config- 
urations. Snapshots of time evolution, phase projection 
and localized sets of the circuits as observed from the 



oscilloscope confirm the existence of different synchro- 
nized states experimentally along with corresponding nu- 
merical results. Further, the transition to different syn- 
chronized states can be numerically quantified from the 
changes in the largest Lyapunov Exponents (LEs), Cor- 
relation of Probability of Recurrence (CPR) and Correla- 
tion Coefficient (CC) of the coupled systems as a function 
of the coupling strength. We also present a detailed linear 
stability analysis and obtain synchronization conditions 
for different synchronized states in the present system. 

The remaining paper is organized as follows: In Sec. HH 
we will describe briefly about the system employed to 
demonstrate different synchronized states and explain 
the circuit realization. In Sec. [Ill] we explain the mutual 
coupling configuration and present the experimental and 
numerical results to confirm various types of synchroniza- 
tions along with the necessary linear stability analysis. 
In Sec. [Vj we describe the subsystem coupling configu- 
ration and demonstrate the occurrence of different syn- 
chronization transitions and, finally, we summarize our 
results with conclusion in Sec. [Vj 



II. DYNAMIC ENVIRONMENT COUPLING 

A. System description 

First we consider a single scalar delay differential equa- 
tion given as 



x = -ax(t)+/3f(x(t-r)) 



(1) 



where the nonlinear function f{x), is represented by a 
piecewise linear function 



f{x) = AF* - Bx. 



Here 



F* — 



—x 
x\ 



x < -x* 

-x* < X < X* 

X > X*. 



(2) 



(3) 



The system parameters are fixed as a = 1.0, j3 = 1.2, 
time-delay r = 6.0, A = 5.2, B = 3.5 and x* = 0.7 
throughout the manuscript. For these chosen parame- 
ter values the single system (pp) exhibits a hyperchaotic 
attractor with three positive LEs [24|, [25j . 



B. Circuit realization 

The design of the electronic circuit which describes the 
dynamics of Eq.(pQ) along with the threshold piecewise 
linear function f(x) is given in Fig. [TJ This circuit con- 
sists of a diode based nonlinear device (ND) with two 
amplification stages {OA\ and OA2), a time-delay unit 
(DELAY) and a low pass first order RqCq filter. Here 
liA74As are engaged as operational amplifiers. V\ and 
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FIG. 2: Circuit block diagram of the coupled time delayed 
feedback oscillator for mutual coupling configuration J4}. 
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V2 are the constant voltage sources of all active devices 
(±12V). Using thses voltage values Vi and V2, one can 
easily adjust the threshold values of the three segment 
piecewise function (Eq. (J3J). By applying the KirchhofPs 
laws to this circuit the state equation can be written as 
RoCo^p- = -U(t) + F[k f (U(t-T d ))], where U(t) is the 
voltage across the capacitor Co, U(t — Td) is the voltage 
across the delay unit, Td = n\jLC is the time-delay, n 
is the number of LC filter units, and F[kf(U(t — Td))] is 
the static characteristic of the ND unit. 

To study the circuit equation, we transform it to the 
dimensionless equation (pQ) by defining the dimensionless 
variables and parameters as x(t) = -^, t' = R t c , 

r = R d c , kf = /3, and t' —> t. The circuit param- 
eters are fixed as i?i = liffi, R2 = R3 = 10KQ, 
R A = 2KQ, R 5 = 3KQ, R 6 = 10 AKVt (trimmer-pot), 
R 7 = IKVL, R 8 = 5KQ (trimmer-pot), (R 9 = R 10 = 
1KQ, Ru = 10KQ, R12 = 20KQ (trimmer-pot), R = 
1.86KQ, C = lOOnF, Li = YlmE (i = 1,2,- •• ,11), 
d = 470nF (i = 1, 2, • • • , 10), n = 10. T d = 0.751ras, 



R0C0 = 0.268ms, and so the time-delay r 
chosen circuit parameter values. 



2.8 for the 



III. MUTUAL COUPLING CONFIGURATION 

We now construct a circuit which consists of a system 
of two identical time-delayed sub circuits with thresh- 
old piecewise linear nonlinearity and are coupled indi- 
rectly through a common environment. Here both the 
circuits and environment are mutually sharing their feed- 
back with each other and the state equation for the cou- 
pled circuit can be written as 



RqCq 



dU^t) 



RqCq 



dt 
dU 2 (t) 



dt 



-a'U^ + ftkfU^t-Td)} 



-a'U 2 (t) + f[k f U 2 (t-T d 



e[f3 2 V(t), 
dV(t) 



(4a) 



(4b) 



RqCq 



dt 



~k'V(t) - -f [P'Mt) + $l/ 2 (t)](4c) 



where U\ (t) and U 2 (t) correspond to the output variables 
of each circuit abd V(t) is the output of the environmen- 



FIG. 3: (Color online) Snap shots of the time evolution of the 
coupled systems (yellow Ui(t) and green t/2 (t) ) with mutual 
coupling configuration J4} obtained from the oscilloscope in- 
dicating the existence of (a) IPS and (b) IS in time-delayed 
electronic circuits. Vertical scale 2.0v/div and horizontal scale 
2.0ms/div. 



tal equation. The schematic diagram for this coupled 
circuit is sketched in Fig. [2j By defining the normalized 



Ul,2(t) 

U s 



variables and parameters as above and x\ :2 (t) 

y(t) = -^ one can obtains the equivalent dimensionless 
equation as follows: 



Xi(t) = 


-axi(t) + fif[xi(t - r)] + eifiiy, 


(5a) 


x 2 (t) = 


-ax 2 (t) + f3f[x 2 (t - t)} + eifcy, 


(5b) 


y = 


s 2 
-ky - —(/3ix 1 + fi 2 x 2 ), 


(5c) 



where, a' = a, /3[ 2 = /3i 5 2, s[ 2 = £1,2 and k' = k. 
Here the two systems x\(t) and x 2 (t) are not directly 
coupled to each other, instead they are coupled through a 
coupling function (y) which has a dynamics modulated by 
the system dynamics. S\ is the strength of the feedback 
to the systems and e 2 is the strength of the feedback 
to the environment (coupling parameters). j3\ and j3 2 
are the nature of feedback from and to the environment, 
respectively, k is the damping parameter and we choose it 
as k = 1. In the absence of feedback from the systems to 
the environment, the strength of the environment decays 
exponentially fast as k > 0. 



A. Experimental and numerical observations 

We have experimentally observed different types of 
synchronization in the coupled electronic circuits which 
are also confirmed using numerical simulations. Depend- 
ing upon the feedback strength /3\ and f3 2 we observe 
two types of synchronization transitions in the coupled 
systems. When f3\ and j3 2 are of the same sign, for exam- 
ple (/?i,/?2) = (1, 1), we observe a transition from non- 
synchronization to IS via IPS and when f3\ and f3 2 are 
of different sign, (/?i,/?2) — (1? — 1)? then we observe a 
transition from non-synchronization to CS via PS as a 
function of the coupling strength. 
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FIG. 4: (Color online) The numerically obtained time series 
for (/?i,/?2) = (1, 1) in mutual coupling configuration (0. (a) 
IPS for e — 0.6 and (b) IS occurs for e — 1.5. 



First, we consider the case of (/?i,/?2) = (1, 1) and 
for simplicity we have chosen the value of the coupling 
strengths E\ = £2 = £• In the absence of the coupling 
(e = 0) both the circuits oscillate independently, and for 
sufficiently large value of the coupling strength the phase 
difference between the circuits is exactly 7r, that is the 
systems exhibit IPS. On increasing the coupling strength 
to further larger values, IS occurrs between the circuits. 

Snapshots of the wave forms of the circuits as seen in 
the oscilloscope are shown in Fig. [3j The existence of IPS 
between the circuits is depicted in Fig. [3fa) where the 
circuits are evolving with phase difference of tt but still 
the amplitudes are uncorrelated. Figure. [3]^b) shows the 
realization of IS between the systems where both phase 
and amplitude are correlated and occur exactly opposite 
to each other. If we numerically solve Eq. (|5]), we obtain 
IPS for e = 0.6 which is depicted in Fig. [4] and this syn- 
chronization can also be confirmed both experimentally 
and numerically using the phase projection plots which 
are shown in Figs. [5fa) and[5fc), respectively. Further 
increase the coupling strength to e = 1.5 the systems 
exhibit IS where the maxima of both systems occur ex- 
actly opposite to each other as depicted in Fig.HJb). The 
corresponding experimental and numerical phase projec- 
tion plots of the systems are given in Figs.[5fb) and[5fd), 
respectively. 

Phase coherence of the systems is further qualitatively 
visualized both experimentally and numerically using the 
framework of localized sets [26|. The basic idea of this 
characterization is that the set of points obtained by sam- 
pling the time-series of the system 1 whenever the maxi- 
mum occurs in system 2 is plotted along with the attrac- 
tor of the system 1 and vice versa. If the coupled systems 
are said to be phase synchronized then the sets are local- 
ized on the attr actor, otherwise they spread over the en- 
tire attractor implying asynchronization. This provides 
an easy and efficient way to detect phase synchronization 
even in non-phase-coherent attractors. 

The experimental observation of the localized sets is 
shown in Fig. [6] In the absence of the coupling [e — 0), 
the sets are distributed over the entire attractor which 
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FIG. 5: (Color online) Experimental (vertical scale 2.0v/div 
and horizontal scale 2.0v/div) and numerical observations of 
phase projections of the coupled systems for mutual coupling 
configuration, (a), (c) IPS occurs for e — 0.6 and (b), (d) IS 
occurs for e — 1.5. 




FIG. 6: (Color online) Experimental realization of the frame- 
work of localized sets, (a), (b) The sets are spread over the 
attractors indicating the absence of phase coherence in the 
absence of the coupling, (c), (d) For a sufficiently large value 
of coupling strength, the sets are localized on the attractors 
which indicates IPS for mutual coupling configuration. Ver- 
tical scale 2.0v/div and horizontal scale 0.5v/div. 



corresponds to the absence of phase coherence. Figures 
[6ja) and[6fb) show that the attractors of the two systems 
along with the sets for the case of nonphase synchroniza- 
tion. The corresponding numerically obtained figures are 
plotted in Figs. [7Ja) and[7Jb) for e = 0. For sufficiently 
large value of the coupling strength, the sets are localized 
on their corresponding attractors which confirm a perfect 
phase locking of the systems [Figs. [6^c) and[6fd)]. The 
corresponding numerical figures are plotted in Figs. [T^c) 
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FIG. 7: (Color online) Numerically obtained corresponding 
figures for the localized sets, (a), (b) The sets are spread over 
the entire attractors in the absence of the coupling (e = 0). 
(c), (d) For the value of e = 0.6, the sets are localized on the 
attractors. 
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FIG. 8: (Color online) (a) CPR (continuous line), CC (dotted 
line) and (b) spectrum of maximal LEs for (/3i,/?2) = (1,1) 
in mutual coupling configuration as a function of coupling 
strength e E (0,2.5). 



and [71(d) for the value of s = 0.6. 

Next, the transition from nonsynchronization to IS via 
IPS can be characterized by the changes in the spectrum 
of maximal LEs. Also the phase coherence is further 
quantified using the index CPR. Complete and inverse 
synchronizations can be quantified using the CC. these 
are given by the expressions, 



CPR = (Pi{t)P 2 {t))/a 1 <j2, 
cc _ ((x^t) - (x 1 (t)))(x 2 (t) 



<*2(*)»> 



(6) 



v /((x 1 (t)-(x 1 (tW)((x 2 (t)-(x 2 (tW 

where ( ) brackets indicate time average. Using Eq. 
(|6j), we calculate the index CPR. Here P\^ means that 




FIG. 9: (Color online) Snap shots of the time evolution of 
both coupled circuits (yellow Ui(t) and green U2(t)) indi- 
cating the existence of (a) PS and (b) IPS in coupled time- 
delayed electronic circuits for mutual coupling configuration 
with (/3i, #2) = (1,-1). Vertical scale 2.0v /div and horizontal 
scale 1.0ms /div. 
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FIG. 10: (Color online) The time series of the coupled systems 
in mutual coupling configuration J5J for (^1,^2) = (1, — 1)- 
(a) PS for e = 0.6 and (b) CS for e = 1.5. 



the mean value has been subtracted, ai^ are the stan- 
dard deviations of P\ and P^ and P(t) is a generalized 
autocorrelation function based on recurrence properties 
|27j . If the phase of the systems are perfectly locked, 
then the probability of recurrence is maximal at the time 
t and CPR ~ 1, otherwise the maxima do not occur si- 
multaneously and hence one can expect a drift in both 
the probability of recurrence resulting in low values of 
CPR. Using Eq. flZ]), we calculate the CC to characterize 
the CS and IS between the systems. If both systems are 
in CS state then, CC « 1 and for IS state CC will be 
« -1. 

In Fig. Ufa) we have plotted the numerically calculated 
CPR (continuous line), CC (dotted line) and in Fig. [51(b) 
we have plotted the ten largest LEs of the coupled sys- 
tems as a function of the coupling strength e E (0, 2.5). In 
the absence of the coupling (e = 0), the index CPR and 
CC are near to zero which confirms that the systems are 
evolving without any synchronization with six positive 
LEs (three for each systems). If the coupling strength 
increases, the index CPR also starts to increase towards 
unit value and for e = 0.56, CPR is oscillating near unity 
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FIG. 11: (Color online) Experimental (vertical scale 2.0v /div 
and horizontal scale 2.0v /div) and numerical phase projection 
plots of the coupled systems with mutual coupling configura- 
tion for (#l,#0 = (1,-1). (a), (c) PS for s = 0.6 and (b), 
(d) IS for e = 1.5. 



(CPR ~ 0.97) whereas the value of CC becomes nega- 
tive which indeed confirms the onset of IPS. Additional 
confirmation comes from the changes in the LEs, where 
the zeroth LE of the coupled system becomes negative 
at e = 0.56 which indicates the existence of IPS. If we 
increase the coupling strength further, the CC of the cou- 
pled systems decreases and reaches the value « —0.99 at 
e = 1.39. At this value of £, except for the three positive 
LEs, all the other positive LEs of the coupled systems 
become negative which indeed confirms the existence of 
IS in the coupled time-delay systems. 

Next, we consider the case with (/^i,/^) = (1,-1) 
(different signs). Now we observe the transition from 
nonsynchronization to CS via PS as a function of the 
coupling strength. In this case, for lower values of cou- 
pling strength the individual circuits evolve indepen- 
dently, while for e = 0.56 the circuits exhibit PS and for 
further larger value of e {e = 1.39) both systems attain 
a CS state. 

The experimental snap shots and numerical plots of 
the time series of the coupled systems are depicted in 
Figs. [9^a) and [TOt a), respectively, for e — 0.6 indicating 
that the circuits are evolving with PS. This is also verified 
(both experimentally and numerically) using the phase 
projection plots [Figs. fTTT a) and fTTT c)], respectively. On 
further in increase the coupling strength to £ = 1.5, the 
circuits exhibit CS with each other which is depicted in 
Figs. [9jb) and fTOT b) (both experimentally and numeri- 
cally, respectively) . The corresponding experimental and 
numerical phase projection plots of the systems are given 
in Figs. HTf b) andQjJd), respectively. 

PS is further confirmed by the frame work of localized 




FIG. 12: (Color online) Experimental realization of the frame- 
work of localized sets in the case of (^1,^2) = (1,-1) in mu- 
tual coupling configuration, (a), (b) The sets are spread over 
the attractors indicating the absence of phase coherence, (c), 
(d) The sets are localized on the attractors which indicates 
phase synchronization. Vertical scale 2.0v/div and horizontal 
scale O.bv/div. 
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FIG. 13: (Color online) Numerically obtained equivalent fig- 
ures for Fig. 1121 (a), (b) The sets are spread over the entire 
attractors for s = and (c,d) the sets are localized on attrac- 
tors which shows PS for e — 0.6. 



sets. The experimental observations of localized sets are 
presented in Fig. [12] The sets are distributed over en- 
tire attractors which confirm no synchronization in the 
absence of the coupling [Figs. fT2fa) andH2Tb)]. The cor- 
responding numerical figures are shown in Figs. fl3T a) 
and fT3fb) for £ = 0. For sufficiently large value of cou- 
pling strength, the sets are localized on the attractors 
which confirm a perfect phase locking between the sys- 
tems [Figs. \T2\c) and[T2Td)]. The equivalent numerical 
figures are given in Figs.fT3Tc) andfl3Td) for the value of 
coupling strength e = 1.5. 

In Fig. EJa) we have plotted the index CPR, CC and 
in Fig. [14] (b) we have plotted ten maximal LEs of the 
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FIG. 14: (Color online) (a) CPR (continuous line), CC (dot- 
ted line) and (b) spectrum of maximal LEs for (/?i,/?2) = 
(1,-1) in mutual coupling configuration as a function of cou- 
pling strength e G (0, 2.5). 




FIG. 15: (Color online) (a), (b) Two parameter diagrams in 
(e\ — £2) plane shows various types of synchronization states 
in mutual coupling configuration for (^1,^2) = (1,1) and 
(^1,^2) = (1, — 1) , respectively. In both figures, the white 
color indicates desynchronized state, red color represents IPS 
and PS regime and green color belongs to IS and CS regime. 
The black points indicate the analytically obtained stability 
condition of the outer regime of the piecewise linear function 
[Eq. 0] and blue (light gray) points indicate the stability 
condition for the middle region of the piecewise linear func- 
tion. 



coupled systems as a function of the coupling strength 
(e G (0,2.5)). In the absence of the coupling, the index 
CPR and CC ~ confirming that the systems are inde- 
pendently oscillating without any synchronization. Once 
the coupling strength increases beyond 0.5 (e = 0.56), 
the index CPR increases and reaches a value close to 
unity (CPR « 0.97) and CC becomes « 0.8, confirming 
the onset of PS but still the amplitudes are uncorrelated 
[Fig. [TIT a)]. PS is also confirmed from the changes in the 
LEs where the zeroth LE of the coupled systems become 
negative at e > 0.5 [Fig.lTSTb)]. Further, for e = 1.39, the 
CC of the systems reach the unit value (~ 0.99) which is 
depicted in Fig. [THa) and for this e value, except for the 
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FIG. 16: Circuit block diagram of the coupled time delayed 
feedback oscillator for subsystem coupling configuration (|T6]) . 



three positive LEs, all the other positive LEs of the cou- 
pled system become negative confirming the existence of 
CS [Fig. Gib)]. 

To identify the global picture of the different synchro- 
nization transition regimes we have plotted the phase di- 
agrams in the two parameter plane of coupling strengths 
(61,62)- We use the index CPR to mark the PS and IPS 
regimes. In addition, CS and IS regimes are marked using 
the CC between the systems. The phase diagrams in the 
(si —82) plane for coupled time-delay systems are shown 
in Fig.H^a) andEJb) for (/3i,/3 2 ) = (1,1) and (1,-1), 
respectively. In Figs. [T5l a) and [T5l b) the red color (dark 
gray) corresponds to IPS and PS regimes and the green 
color (light gray) regime indicates IS and CS states. Also 
the white color corresponds to desynchronized state. 



B. Linear stability analysis 

In this section, we give a detailed stability analysis for 
the synchronized states of the two dynamically coupled 
time-delay systems. To find the stability of the synchro- 
nization, we apply a infinitesimal perturbation to the sys- 
tem ([5]) . Let £1 , £2 and z represent the perturbations then 
Eq. (J5j) can be rewritten as 

£1 = g'(xi,xi T )€i+6i/3iz, (8a) 

£2 = g'(x2iX2T)&+6 1 f3 2 Z, (8b) 

hti+fafr]. (8c) 



-kz- S 4\ 



Here g'(x,x T ) is the derivative of the Eq. (J2j) and 
for the two outer regime of Eq. ([3]) it can be written 
as g'(x,x T ) = —(a + /3B) and for the middle region 
g'(x, x T ) = -a + f3(A - B). 

Equation JHJ) is difficult to solve with the present form 
and so, we consider the special case of complete synchro- 
nization state (x\ — X2). Now Eq. (jSJ) can be simplified 
by denning 



So Eq. (|8j) becomes, 

z = -kz- y 7 f . 



(9) 

(10a) 
(10b) 



Here A = g'(x, x T ). The synchronized state is said to be 
stable if all the LEs obtained from Eq. ([TO]) are negative. 
Eq. [10] can be further simplified by choosing fa 1 + /3f = 2 
and eliminating z from it so that we get an equation 



£ = (A - fc)£ + (Afc - £i£ 2 )£o. 
Consider a solution of the form 



Co = Ae r 



so that we get 



(A - k) ± y/(X - fc) 2 - 4(£i£ 2 - Afc) 



(11) 



(12) 



(13) 



The synchronized state is stable if Re[m] is negative for 
both the solutions. This provides the following criteria 
for the stable synchronized state: 

1. If (A — fc) 2 < A(e\E2 — Afc), then m is complex and the 
stability condition is fc > A. 

2. If (A — fc) 2 > 4(£i£2 — Afc), then m is real and the 
stability condition becomes fc > A and 8162 > Afc. Here 
the transition to stable synchronization is given by the 
threshold values of parameters satisfying the condition 



£i > 



(A + fc) 2 

4e 2 



(14) 



From Eq. ([H]) one can obtain the threshold value of the 
coupling strength for the synchronized states. For in- 
stance we choose A = — a+faA — B) corresponding to the 
middle regime (\x\ < x*) of Eq. ^ and A = -(a + f3B) 
corresponding to the outer regime (\x\ > x*) and the 
feedback (fa, fa) = (1? 1) we obtain a transition from no 
synchronization to IS state which is plotted along with 
the numerically obtained synchronized region. The black 
points in Fig. fT5fa) correspond to the condition for the 
stable synchronization state for the outer regime and the 
blue points indicate the stability condition of the middle 
regime of the piecewise linear function. This shows that 
the numerically obtained IS regime occurs in between the 
analytically obtained regime. Similar transition curves 
are also observed for transition from no synchronization 
to CS in the case of (fa, fa) = (1,-1) [Fig. H^b)]. 



IV. SUBSYSTEM COUPLING 
CONFIGURATION 

Further, we have also considered a second form of 
coupling called subsystem coupling configuration, where 
both the circuits share their feedback with the environ- 
ment, while only one of the circuits is receiving the feed- 
back from the environment. The schematic circuit block 
diagram for this coupling configuration is sketched in 









FIG. 17: (Color online) Snap shots of the time evolution of 
both coupled circuits (yellow Ui(t) and green Uztf)) indicat- 
ing the existence of (a) IPS and (b) IS in coupled time-delayed 
electronic circuits for subsystem coupling configuration JTSJ) 
with (fa, fa) = (1,1). Vertical scale 5.0v/div and horizontal 
scale 2.0ms/div. 
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FIG. 18: (Color online) The time series of the coupled sys- 
tems in subsystem coupling configuration (|16|) with (fa, fa) — 
(1, 1). (a) IPS for e = 1.2 and (b) IS for e = 2.0. 



Fig. [TBI and the state equations can be given as follows: 

dU^t) 



RoCo 
RoCo 



dt 

dU 2 (t) 

dt 



e[p f 2 V(t), 



RqCq 



dV(t) 
dt 



= -a'C/ 1 (t) + /[fc / [/ 1 (i-r d )], (15a) 

= -a'U 2 (t) + f[k f U 2 (t-T d )} + 

(15b) 
= -k'V(t)- £ 4iP[Ui(t)+/3 2 U 2 (tlpc) 



and the corresponding dimensionless equations for this 
configuration can be defined as 

x t (t) = -azi(t)+/3/[ffi(£-T)], (16a) 

x 2 (t) = -ax 2 (t)+f3f[x 2 (t-T)]+s 1 fay, (16b) 



V = ~ky- — (fax! +fax 2 ). 



(16c) 



Here the circuits of U±(t) and U 2 (t) (x\(t) and x 2 (t)) 
are sharing their feedback with the environment but only 
the circuit of U 2 (t) is receiving feedback from the en- 
vironment. The parameters are fixed as in the previ- 
ous case. First we consider the case (fa, fa) = (1,1) 
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FIG. 19: (Color online) Experimental (vertical scale 2.0v/div 
and horizontal scale 2.0v /div) and numerical observation of 
phase projection plots for the subsystem coupling configura- 
tion (32} with (0i, #2 ) = (1, !)• (a), (c) IPS for e = 1.2 and 
(b), (d) IS for £ = 2.0. 




FIG. 20: (Color online) Experimental realization of the frame- 
work of localized sets in subsystem coupling configuration 
with /3i,/?2 = (1,1)- (a), (b) The sets are spread over the 
attractors indicating the absence of phase coherence, (c), (d) 
The sets are localized on the attractors which indicates the 
occurrence of IPS. Vertical scale 2.0v /div and horizontal scale 
0.5v/div. 



and we obtain the same kind of transition from IPS (for 
e > 1.0) to IS (for e > 1.9). The experimental and 
numerical snap shots of the time evolution of the cou- 
pled systems show IPS for the value of coupling strength 
£ = 1.2 which is evident from Figs. fTTTa) andlTSTa), re- 
spectively. This can also be confirmed both experimen- 
tally and numerically using the phase projection plots 
[Figs. fT9Ta) and fl9Tc)]. If we increase the coupling to 
further larger values, at e = 2.0 both circuits attain IS as 
depicted in Figs.fTTTb) andfTBTb) (experimentally and nu- 
merically), respectively. The corresponding experimental 
and numerical phase projection plots are presented in 
Figs. [T9Tb) and[T9Td). respectively. 

The localized sets plots are again presented to con- 
firm the existence of IPS. The experimentally obtained 
attractors along with the sets are shown in Figs. l20Ta) 
and l20f b) where the sets are distributed over entire at- 
tractors in the absence of any coupling and the corre- 
sponding numerical figures are depicted in Figs. l2TT a) 
and l2Jlb) for £ = 0. For sufficiently larger values of the 
coupling strength, one can observe that the sets are lo- 
calized on the attractors confirming the perfect locking 
of the phase of the systems [Figs. l20t c) and l20T d)]. The 
equivalent numerical figures are plotted in Figs.EjJc) and 
Eld) for £ = 1.2. 

The transition from IPS to IS can again confirmed by 
plotting the index CPR, CC and the changes in the max- 
imal LEs as a function of the coupling strength. In the 
absence of the coupling the systems are evolving with- 
out any synchronization (with six positive LEs) and so 
CPR and CC w 0. For e = 1.02, the value of index 
CPR increases towards unity and reaches the value CPR 
~ 0.97, whereas CC becomes negative which confirms 




1 x 2 (t) 1 



FIG. 21: (Color online) Numerically obtained equivalent fig- 
ures for Fig. [20] (a), (b) The sets are spread over the entire 
attractors for e — and (c), (d) the sets are localized in one 
part of the attractor which shows IPS for e = 1.2. 



the onset of IPS state [Fig. l22T a)]. This transition is also 
confirmed from the changes in the LEs where the zeroth 
LE of the coupled system becomes negative for this value 
of £ indicating the existence of IPS which is evident from 
Fig. I22fb). On further increase of the coupling strength 
to £ = 1.91 the systems exhibit IS transition where the 
CC becomes w -0.99 [Fig. 122(a)] and from Fig. [22(b), we 
notice that except for the three positive LEs, all the other 
positive LEs become negative confirming the existence of 
IS in the coupled time-delay system. 

Finally, we consider the case (/?i,/?2) — (1? — !)• We 
obtain a transition from no synchronization to CS via 
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FIG. 22: (Color online) (a) CPR (continuous line), CC (dot- 
ted line) and (b) spectrum of maximal LEs for (/?i, /fe) = (1, 1) 
in subsystem coupling configuration as a function of coupling 
strength e € (0,2.5). 
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FIG. 23: (Color online) Snap shots of the time evolution of 
both coupled circuits (yellow U\(t) and green ^(t)) indicat- 
ing the existence of (a) PS and (b) CS in coupled time-delayed 
electronic circuits for subsystem coupling configuration (|15|) . 
Vertical scale 5.0v/div and horizontal scale 1.0ms /div. 



PS in the coupled time-delayed electronic circuit. The 
experimental snap shots of the time evolution of both 
circuits are displayed in Fig. l23T a) indicating that the 
systems exhibit PS and the numerically obtained figure 
is depicted in Fig. IMTa) for e = 1.2. The corresponding 
experimental and numerical phase projection plots are 
displayed in Figs. I25l a) and l25l c), respectively. CS is 
observed between the circuits for further larger value of 
the coupling strength and the experimental wave forms 
for suitable e is shown in Fig.l23Tb) exhibiting CS. Figure 
B^Tb) shows the numerically obtained time traces of the 
two systems displaying CS and Figs. l25f b) and l25f d) 
show the corresponding phase projections (obtained by 
experimental and numerical simulations, respectively) of 
the systems for e = 2.0. 

PS between the circuits can be once again visualized 
by plotting the localized sets. Figure [26] shows the ex- 
perimentally observed attractors along with the sets. In 
Figs. l26l a) and l26f b) the sets are distributed over the 
entire attractor for lower values of coupling strength due 
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FIG. 24: (Color online) Numerically obtained time series of 
the coupled systems in subsystem coupling configuration (|16p 
with (Pi,p 2 ) = (1,-1). (a) PS for e = 1.2 and (b) CS for 
£ = 2.0. 
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FIG. 25: (Color online) Experimental (vertical scale 2.0v/div 
and horizontal scale 2.0v /div) and numerical phase projection 
plots of the coupled systems with subsystem coupling config- 
uration ([16} with (#l,/3 2 ) = (1, -!)• (a), (c) PS for e = 1.2 
and (b), (d) IS for e = 2.0. 



to the absence of PS. The corresponding numerically ob- 
tained figures are plotted in Figs. l27Ta) and [27Tb) for 
e = 0. If we increase the coupling strength to a suffi- 
ciently large value, the sets are localized on the attractor 
as depicted in Figs. l26T c) and l26T d) confirming the phase 
locking of both the systems. The corresponding numeri- 
cally obtained figures are plotted in Figs. l27f c) and l27T d) 
for the value of e = 1.2. 

In the absence of the coupling the systems are evolving 
freely without any synchronization with six positive LEs 
(three for each systems) and so CPR and CC ~ 0. Be- 
yond e > 1.0 (e = 1.02), CPR becomes « 0.97 which con- 
firms the existence of PS which is depicted in Fig.,l28ta). 
At this value of £, the zeroth LE of the coupled sys- 
tem becomes negative which is shown in Fig. l28T b). If 
we increase the coupling strength further, the CC of the 
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FIG. 26: (Color online) Experimental realization of the lo- 
calized sets for the case of (/?i,/?2) = (1,-1) in subsystem 
coupling configuration, (a), (b) The sets are spread over the 
attractors indicating the absence of phase coherence, (c), (d) 
The sets are localized on the attractors which indicates phase 
synchronization. Vertical scale 2.0v/div and horizontal scale 
0.5v/div. 
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FIG. 27: (Color online) Numerically obtained equivalent lo- 
calized sets figures for Fig. [26l (a), (b) The sets are spread 
over the entire attractors for e — and (c), (d) the sets are 
localized in one part of the attractor which confirms IPS for 
s = 1.2. 



systems increases and reaches the unit value at e = 1.91 
which is evident from Fig. l28T a) and for this value of cou- 
pling strength all the positive LEs (except three positive 
LEs) become negative confirming the existence of CS in 
the coupled time-delay systems [Fig. ITSTb)]. 

The phase diagrams in the parameter plane of coupling 
strengths are plotted to identify the global picture of the 
regimes of different types of synchronization states. The 
phase diagrams in the (s\ — £2) plane for coupled time- 
delay systems (Eq.fT6|) are shown in Fig. [29t a) and l29f b) 
for (/?i,/?2) = (1, 1) and (1,-1), respectively. In these 




FIG. 28: (Color online) (a) CPR (continuous line), CC (dot- 
ted line) and (b) spectrum of maximal LEs for subsystem 
coupling configuration with (^1,^2) = (1,-1) in subsys- 
tem coupling configuration as a function of coupling strength 
e e (0,2.5). 




FIG. 29: (Color online) (a), (b) Two parameter diagrams in 
(si — £2) plane shows various types of synchronization states 
in subsystem coupling configuration for (^1,^2) = (1, 1) and 
(^1,^2) = (1 , — 1) , respectively. In both figures, the white 
color indicates desynchronized state, red color represents IPS 
and PS regime and green color belongs to IS and CS regime. 
The black points indicate the analytically obtained stability 
condition of the outer regime of the piecewise linear function 
(Eq. Q) and blue (light gray) points indicate the stability 
condition for the middle region of the piecewise linear func- 
tion. 



figures, red color (dark gray) corresponds to IPS and PS 
regime, respectively, and green color (light gray) indi- 
cate IS and CS states, respectively. Also the white color 
corresponds to desynchronized state. 

We have also carried out a linear stability analysis for 
the above coupling configuration as in Sec. IIIIBI and it 
gives a condition for the stability of the synchronized 



states as e\ > 



(A+fe) 2 
2e 2 



From this relation one can obtain 



the threshold value of the coupling strength for different 
synchronized states. For example, we choose A values as 
in Sec. IIIII and the feedback (/?i,/?2) = (1? 1)? we obtain 
a transition from no synchronization to IS state which is 
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plotted along with the numerically obtained synchronized 
region. The black points in Fig. l29T a) correspond to 
the outer regime of the piecewise linear function ((3]) and 
the blue points indicate the stability state of the middle 
regime. Similar to Figjl5j the numerically obtained syn- 
chronization regime occurs in between the analytically 
obtained regimes. Figure l29f b) shows the transition 
curves for the case of (/?i,/?2) = (1, — 1) which indicates 
the transition from no synchronization to CS in subsys- 
tem coupling configuration. 



V. CONCLUSION 

In this paper, we have experimentally demonstrated 
the occurrence of various types of synchronization in 
coupled piece-wise linear time-delayed electronic circuits 
with threshold nonlinear function where the circuits are 
coupled indirectly through a common dynamic environ- 
ment. We have carried out these studies in two different 
coupling configurations, namely mutual and subsystem 
coupling configurations. In both configurations, depend- 
ing upon the strength of the coupling and the nature of 
the feedback, we observe different types of synchroniza- 
tion transitions in the coupled circuits such as transi- 
tion from IPS to IS and from PS to CS in hyperchaotic 



regimes. Snapshots of the time evolution, phase projec- 
tion and localized sets plots of the circuits observed from 
the oscilloscope confirm the various synchronization phe- 
nomenon experimentally. The corresponding numerical 
simulations are also presented in detail. Further, the 
transition to different synchronization states can be ver- 
ified from the changes in the maximal LEs, index CPR 
and CC of the coupled systems as a function of the cou- 
pling strength. Also we have presented a detailed linear 
stability analysis to obtain synchronization conditions for 
different synchronized states. 

Note: After the completion of this work, we have 
come across the recent arXiv paper by T. Banerjee and 
D. Biswas J28j, which deals with a similar problem on 
synchronization in hyperchaotic time-delayed electronic 
oscillators coupled indirectly via a common environment 
using a tank nonlinearity under bidirectional coupling, 
confirming the occurrence of in-phase and inverse-phase 
synchronization phenomena. 
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